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A  Statistical  Solution  to  Shallow  Water  Propagation  that  Includes  Chaos 

N00014-90-J>13d5 

Martin  A.  Mazur  and  Diana  F.  McCammon 
The  Pennsylvania  State  University 
Applied  Research  Laboratory 
P.  O.  Box  30 
State  College,  PA  16804 

Tel:  (814)  863-4782  (Mazur)  (814)  863-4159  (McCammon) 
email:  internet:  mam@arlvax.arLpsu.edu 

Objective:  Ray  tracing  in  range  dependent  environments  generally  leads  to  chaotic  ray  paths 
[1].  However,  Collins  and  Kuperman  [2,3]  have  recently  shown  that  ray  chaos  can  be  overcome 
by  recasting  ray  tracing  in  terms  of  Fermat’s  principle  of  minimum  propagation  time.  The 
objective  of  this  research  is  to  expand  this  idea  by  flying  a  variety  of  techniques  from 
optimization  theory  and  the  theory  of  two-point  boundary  value  problems  to  the  calculation  of 
dgenrays  in  range  dependent  environments. 

Background:  In  classical  ray  tracing,  eigenrays  between  a  source  and  receiver  are  determined 
by  an  initial  value  or  ‘shooting’  tq>proach.  The  launch  angles  of  rays  from  a  source  point  are 
varied  until  the  rays  intersect  the  receiver  end  point.  In  range  dependent  environments,  the  ray 
paths  are  generally  chaotic  [1].  This  puts  a  fundamental  limit  on  the  accuracy  of  ray  tracing, 
evoi  with  perfect  knowledge  of  environmental  conditions.  Because  of  the  extreme  sensitivity 
of  chaotic  ray  paths  to  initial  conditions,  the  calculation  of  eigenrays  in  range  depoident 
environments  quickly  becomes  impossible.  Recently,  Collins  and  Kuperman  [2,3]  have  shown 
that  ray  chaos  can  be  overcome  by  recasting  the  problem  in  terms  of  Fermat’s  principle  of 
minimum  pnqngaticm  time.  The  pi^lem  then  becomes  amenable  to  techniques  of  optimization 
theory  and  the  theory  of  two-point  boundary  value  problems.  In  [3],  the  travel  time  integral  was 
minimized  using  an  optimization  technique  c^ed  simulated  annealing.  A  relaxation  technique 
for  the  solution  of  two-point  bounda^  value  problems,  the  ‘bending’  method,  was  also 
investigated.  Simulated  annealing  produces  stable  eigenrays,  but  is  computationally  inefficiait. 
The  rdaxation  technique  is  very  fast  but  extremely  sensitive  to  the  quality  of  the  initial  guess 
and  can  fail  to  convoge. 

Cnrrent  Researdi:  AccompUshments  and  Goals.  In  our  research  we  are  investigating  using 
the  Rayteigh-Ritz  method  to  minimize  the  travel  time  int^ral.  In  this  technique,  a  candidate 
eigenray  is  rqnesented  as  a  finite  linear  combination  from  among  a  complete  set  of  basis 
functions.  The  coefficients  of  these  functions  are  adjusted  until  the  travel  time  is  minimized. 
The  number  of  coefficients  is  then  increased  until  the  desired  accuracy  is  achieved.  Using  this 
technique,  we  have  calculated  stable  eigenrays  in  a  range  dependent  sound  speed  profile  that  is 
notable  for  producing  extremely  chaotic  ray  paths  when  the  shooting  method  is  used.  Figure 
1  dtows  some  rays  emanating  ftom  a  source  point  as  calculated  using  the  shooting  method.  The 
launch  angles  of  these  rays  are  all  within  lO*^  degree  of  the  launch  angle  of  the  eigenray 
connecting  (0,0.2)  and  (100,  0.2)  calculated  using  simulated  annealing  [3].  The  eigenray 
calculated  using  the  Rayleigh-Ritz  method  is  indistinguishable  from  this  eigenray  and  is  shown 
in  bold.  The  sound  spe^  profile  is  c(r,z)  =  Cq  (1  +  a  cos  K,  r  cos  K^)  where  Cq  =  1.5  km/s, 
a  -  0.02,  Kj  »  2ir/0.5  km**  and  K,  =  2t/20  km  *. 
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Figure  1.  Rayleigh-Ritz  Eigenray  and  Chaotic  Shooting  Rays 

We  have  used  both  sine  functions  and  ‘hat’  functions  as  basis  functions.  The  sine  function  basis 
set  is  particularly  useful  for  the  following  reasons:  (1)  The  resulting  representation  for  an 
eigenray  is  economical.  In  our  example,  the  eigenray  is  represented  well  by  16  coefficients. 
This  compares  with  having  to  optimize  over  hundreds  of  grid  points  in  the  simulated  annealing 
and  relaxation  techniques.  (2)  In  environments  where  the  rays  oscillate,  the  coefficients  found 
in  the  calculation  of  an  eigenray  can  be  useful  in  formulating  initial  guesses  for  the  coefficients 
of  longer  range  rays. 

Our  immediate  goal  is  to  submit  for  publication  in  JASA  a  paper  detailing  our  results  to  date. 
Mr.  Mazur  will  continue  the  research  started  under  this  grant  to  complete  his  doctoral  thesis. 
This  research  will  include  investigations  of  other  optimization  techniques  and  comparisons  of 
their  usefulness  in  terms  of  ease  of  implementation,  speed,  and  adaptability  to  environments  with 
boundary  interactions. 

References 

[1]  D.  R.  Palmer.  M.  G.  Brown.  F.  D.  Tappert.  H.  F.  Bezdek,  "Classical  Chaos  in 
Nonseoarable  Wave  Propagation  Problems".  Geoph.  Res.  Let.,  v.  15  No.  6.  569-572 
(1988). 

[2]  W.  A.  Kuperman  and  M.  D,  Collins.  "Finding  eigenrays  by  optimization  with 
application  to  tomography  and  overcoming  chaos".  J.  .\coust.  Soc.  Am  93  2425  (1993). 

[3]  M.  D.  Collins  and  W.  A.  Kuperman.  "Overcoming  ray  chaos",  to  appear  in  June  1994 
J.  Acoust.  Soc.  .\m. 
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ConfanNKCS  Attended  and  Contributed  Talks:  Martin  A.  Mazur  attended  the  124th  Meeting 
of  the  Acoustical  Society  of  America  held  in  New  Orleans,  Louisiana  from  October  31  to 
November  4,  1992.  He  delivered  a  talk  related  to  this  research  effort  entitled  Does  ray  chaos 
in  range  dependent  environments  disappear  when  higher-order  approximations  are  made? 


Graduate  Students:  Martin  A.  Mazur  is  a  graduate  student  supported  by  this  research  effort. 
He  is  a  US  citizen.  His  expected  date  of  graduation  if  June,  1995.  The  title  of  his  Ph.D.  thesis 
is  Novel  Approaches  to  Mitigating  the  Effects  of  Chaos  in  the  Calculation  of  Acoustic  Eigenrays 
in  Range  Dependent  Environments.  The  thesis  will  be  theoretical  with  significant  numeric^ 
portions.  The  objective  of  the  thesis  is  to  develop  and  compare  a  variety  of  techniques  which 
can  be  used  to  mitigate  or  overcome  ray  chaos  when  calculating  eigenrays  in  range  dependent 
environments.  The  theoretical  underpinnings  of  these  techniques  will  also  be  developed. 
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High  Frequency  Surface  Scattering  Statistics  (N00014-90-J-1365) 

Richard  Lee  Culver 
The  Pennsylvania  State  University 
Graduate  Program  in  Acoustics 
State  College.  PA  16804 
(814)  867-0150 
Email:  rlc@abe.arl.psu.edu 


FY93  Accomplishmoits 

•  Completed  development  and  testing  of  a  five  axis  transducer  positioning  system  for  use 
in  the  new  29’x23’xl8’  acoustic  measurement  water  tank  (pictured  above). 

•  Designed  and  built  a  very  high  speed  (1  ms)  transmitAeceive  (T/R)  switch  for  monos¬ 
tatic  acoustic  measurements. 

•  Non-Rayleigh  scattering  encountered  using  a  single  frequency  sinusoidal  surface  when 
returns  from  several  aspect  angles  are  combined. 
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Papers,  Talks  and  Conferences 


•  P.D.  Neumann  and  R.L.  Culver,  Statistical  characteristics  of  bistatic  surface  scatter. 
presentation  made  at  the  Acoustical  Society  of  America  meeting,  Ottawa,  Canada,  17-21 
May  1993. 

•  Symposium  on  Coastal  Oceanography  and  Litmral  Warfare,  San  Diego,  CA,  2-5  August 
1993. 

•  R.L.  Culver  and  P.D.  Neumann,  High  Frequency  Bistadc  Ocean  Surface  Scattering, 
invited  paper.  Acoustical  Society  of  America  semi-annual  meeting,  E>enver,  Colorado, 

4-8  October  1993. 

•  K.M.  Becker,  Construction  and  calibration  of  an  underwater  acoustic  transducer, 
invited  article.  Scientific  American,  mailed  February  1994. 

•  R.L.  Culver  and  K.M.  Becker,  A  computer-controlled  5  axis  transducer  positioning  sys¬ 
tem,  under  preparation. 

Graduate  Students  Supported 

Kyle  M.  Becker  (US  citizen)  Expected  graduation  date:  12/94  (MS  Acoustics) 

Experimental/theoretical  thesis:  The  relationship  between  surface  characteristics 

and  scattering  strength  statistics 

Thesis  objective:  The  objective  of  Kyle’s  work  is  to  determine  the  physical  mechanisms  or 
surface  properties  which  cause  ocean  surface  scattering  strength  distribution  to  deviate  from 
expmientiaL  As  part  of  his  research,  he  designed,  built  and  tested  the  transducer  positioning 
system.  He  is  currently  designing  random  but  statistically  stationary  surfaces,  measuring  the 
statistics  of  monostatic  scattering  from  the  surfaces,  and  relating  the  elevation  statistics  of  the 
surfaces  to  the  scattering  statistics. 


Todd  A.  Mathias  (US  citizen)  Expected  graduation  date:  12/95  (MS  Acoustics) 

Experimental/theoretical  thesis  A  Stutt'  of  Bistatic  Surface  Scattering  Statistics 

Thesis  objective:  The  objective  of  Todd’s  work  is  to  investigate  high  frequency  bistatic  sur¬ 
face  scattering  strength.  As  part  of  his  research,  Todd  is  developing  the  computer  interface  to 
the  transducer  positioning  system.  During  the  year  ahead,  Todd  will  extend  Kyle’s  work  to 
bistadc  geometries  (tran^tter  and  receiver  not  at  the  same  location)  and  utilize  a  higher  fre¬ 
quency  transducer  in  order  to  investigate  the  effects  of  scale. 
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Transitions,  Application  of  this  Research 


Because  I  also  direct  the  6.2  Torpedo  Environments  project  (RO35H01)  which  provides 
high  frequency  environmental  acoustic  information  to  the  Advanced  Undersea  Weapons  G&C 
Block  (OT3A),  I  am  able  to  focus  my  6.1  Fundamental  Initiatives  research  on  an  issue  which 
is  critically  important  to  shallow  water  undersea  weapon  development,  non-Rayleigh  boundary 
scattering.  On  the  other  hand,  I  am  also  in  a  position  to  integrate  results  of  the  6.1  research 
inU)  a  6.2  program  for  application  and  evaluation,  and  when  warranted,  transition  to  6.3 
undersea  we^ns  or  MC^  programs. 

My  Fundamental  Initiatives-sponsored  work  is  only  3  years  old,  and  during  the  last 
three  years  we  have  interrupted  our  measurement  program  to  develop  important  tools.  A 
linear  actuator-based  transducer  positioning  system  capable  of  translation  in  3  dimensions  and 
rotation  in  azimuth  and  elevation  was  design^  and  built  because  meaningful  false  alarm 
statistics  require  10^  to  10^  independent  measurements  for  each  geometry.  A  transmit/receive 
(T/R)  switch  (which  blocks  out  the  power  amplifier  during  receive)  capable  of  1  ms  switching 
speeds  was  designed  and  built  because  such  a  device  is  not  commerci^y  available.  A  40  dB 
gain  in  SNR  was  achieved.  Finally,  the  work  has  been  moved  to  a  29*x23'xl8’  water  tank 
located  in  the  new  ARL  building,  completed  only  this  past  faU.  While  the  old  tank  was  ade¬ 
quate,  the  new  tank  provides  a  greater  opportunity  for  the  work  to  expand. 

The  significant  outcome  of  the  work  thus  far  has  been  the  finding  that  non-Rayleigh 
statistics  ate  obtained  when  scattering  strengths  obtained  from  different  aspects  of  a  strongly 
polarized  surface  ate  combined.  This  has  led  to  speculation  that  non-Rayleigh  scattering  from 
the  ocean  boundaries  might  be  caused  by  patches  of  strong  (or  coherent)  scattering  intermin¬ 
gled  with  areas  which  scatter  more  wealdy.  Such  surfaces  have  been  dubbed  corked- 
coherent  by  DJ^.  McCammon,  who  was  a  co-PI  when  this  work  was  began  and  continues  to 
be  invrdved 

The  current  effort  is  to  fabricate  a  series  of  new  surfaces  and  evaluate  the  confused- 
coherent  hypothesis.  Dr.  Steve  Stanic  of  NRL/Stennis  has  agreed  to  provide  information 
about  bottom  relief  obtained  from  his  experiments,  which  will  guide  selection  of  the  surface 
characteristics.  A  multi-axis  milling  machine  at  ARL  can  form  surfaces  with  specified  eleva¬ 
tion  qiectra.  Measurements  will  be  made  with  both  monostatic  and  bistatic  geometries 

This  work  could  one  day  be  applied  to  undersea  weapon  development  in  the  following 
way.  It  is  known  that  boundt^  scattering  can  mask  target  echoes  in  shallow  water.  Our 
w<nk  to  characterize  non-Rayleigh  surface  scattering  can  help  us  understand  the  statistics  of 
shallow  water  boundary  reverberation.  An  effon  under  die  6.2  Torpedo  Environments  Project 
will  atteiiq)t  to  distinguish  target  echoes  from  boundary  reverberation  based  upon  their 
different  statistics. 
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A  Higher  Order  Stochastic  Model  for  Surface  Backscatter  in  High  Sea  States 

N00014-90-J-136S 

Kenneth  E.  Gilbert  and  Lucy  J.  Kulbago 
Applied  Research  Laboratory  and  the  Graduate  Program  in  Acoustics 
The  Pennsylvania  State  University 
P.  O.  Box  30 
State  College,  PA  16804 
Telephone:  (814)  863-8291 

Objective:  To  incorporate  into  an  existing  plane-wave  surface  scattering  model  the 
improvements  necessary  to  model  backscatter  in  high  sea  states  and  shallow  water  where  the 
insonifying  field  departs  significantly  from  a  plane  wave.  Work  to  date  has  dealt  with  high  sea 
states  since  experimental  data  is  available. 

Background:  In  the  previous  year  a  stochastic  oceanic  bubble  layer  model  was  developed  in 
which  the  statistics  of  the  microbubble  distribution  are  controlled  by  turbulence  in  the  inertiai 
subrange  (Kolmogorov  turbulence).  In  FY93  surface  backscatter  was  computed  using  the 
stochastic  bubble  layer  model  and  a  first-order  scattering  method  (Plane-Wave  Bom 
Approximation-PWBA)^  as  well  as  a  higher-order  scattering  method  (Distorted- Wave  Bom 
Appioximation"DWBA)^.  The  results  indicate  that,  in  high  sea  states,  turbulent  transport  of 
microbubbles  has  two  major  effects  on  the  backscatter  from  the  bubble  layer.  First,  the  chaotic 
nature  of  turbulence  leads  to  a  sound  speed  that  varies  stochastically  in  space  and  time  but  that 
has  a  characteristic  wave-number  spectrum.  As  a  consequence,  the  backscatter  from  the  layer 
has  a  characteristic  dependence  on  the  frequency.  Second,  the  bubble  layer,  although  stochastic, 
is  a  region  where  the  sound  speed  is  on  average  less  than  in  bubble-free  water.  In  high  sea 
states  the  average  sound  speed  reduction  is  sufficient  to  cause  significant  near-surface  upward 
refraction  for  low  grazing  angles.  In  such  situations,  the  insonifying  field  is  not  well  represented 
by  a  plane-wave.  Therefore  we  have  considered  a  higher-order  scattering  model  that  takes  into 
account  both  the  average  and  stochastic  properties  of  the  microbubble  distribution.  It  is  shown 
that  with  measured  oceanogriq)hic  inputs  (i.e.,  no  adjustable  parameters)  the  model  gives  a  good 
account  of  the  measured  backscatter  data  of  Ogden  and  Erskin^. 

Approach:  To  account  for  a  non-planar  incident  wave,  we  use  the  so-called  "Distorted-Wave 
Bom  Approximation"  (DWBA).  With  the  DWBA  the  plane-wave  insonifying  field  is  replaced 
by  the  solution  for  a  wave  propagating  in  the  average  sound-speed  profile  in  the  bubble  layer. 
In  Figure  1,  we  compare  the  effective  insonifying  fields  for  plane  waves  (PWBA)  and  distorted 
waves  (DWBA)  for  conditions  typical  of  a  high  sea  state  (WS=  30  kt).  It  is  evident  from 
Figure  I  that  for  low  grazing  angles  and  high  sea  states  there  are  significant  differences  between 
the  plane-wave  solutions  and  the  corresponding  distorted-wave  solutions.  The  effect  of  the 
distortion  of  the  plane  wave  is  to  enhance  the  effective  near-surface  field  at  low  frequencies  (f< 
300  Hz)  and  to  suppress  it  at  higher  frequencies  (f  ^  300  Hz). 

Results:  Figure  2  compares  the  predicted  scattering  strength-versus-frequency  with  the  Ogden- 
Erskin  empirical  fits  to  Critical  Sea  Test  data  for  several  grazing  angles.  Two  distorted  wave 
calculaticms  are  shown:  An  exact  numerical  solution  that  shows  resonance  effects  and  a  WKB 
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frequencies,  this  modification  of  the  plane-wave  result  is  in  qualitative  agreement  with  the 
measured  backscatter.  Current  research  is  direcred  toward  a  better  estimate  for  the  average 
insoni^fing  field  in  order  to  improve  the  agreement  with  measured  surface  backscatter  for  high 
sea  states  and  low  grazing  angles. 

References 


[1]  Kenneth  E.  Gilbert,  A  stochastic  model  for  scattering  from  the  near-surface  bubble  layer, 
J.  Acoust.  Soc.  Am.  94,  3325-2352  (1993). 

[2]  Lucy  J.  Kulbago  and  Kenneth  E.  Gilbert,  "Surface  scattering:  The  role  of  turbulent 
mixing  of  microbubbles,  ”  to  be  published  in  a  refereed  book  of  the  Proceedings  of  the 
International  Conference  on  Theoretical  and  Computational  Acoustics,  Mystic 
Connecticut,  July  5-9,  1(>93. 

[3]  P.  M.  Ogden  and  F.  T.  Erskin.  Naval  Research  Laboratory  Rpt.  NRLyFR/5 160-92-9377, 
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plane  waves  (PWBAl  and  for  distorted  waves  iDWBA). 
Relative  to  the  plane-wave  result,  distortion  of  the  incident 
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Figure  2.  Scattering  strength  versus  frequency  at  several 
flxed  grazing  angles.  The  solid  line  is  a  distorted-wave 
Bom  approximation  calculation  (DWBA)  using  an  exact 
numerical  solution  for  the  distorted  incident  wave.  The 
long  dashed  line  is  a  DWBA  calculation  using  a  WKB 
approximation  for  the  distorted  incident  wave.  The  short 
dashed  line  is  a  plane-wave  Bom  approximation  calculation 
( PWB  A).  The  diamonds  are  the  Ogden-Erskin  empirical  fit 
to  the  Critical  Sea  Test  data. 
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Activity  Report  for  FY93  (1  Oct  1992  -  30  Sept  1993) 


1.  Papers  in  refereed  publications 

Kenneth  E.  Gilbert,  A  stochastic  model  for  scanering  from  the  near-surface  oceanic  bubble 
layer,  J.  Acoust.  Soc.  Am.  94,  3325-2352  (1993). 

Kenneth  E.  Gilben  and  Xiao  Di,  A  fast  Green's  function  method  for  outdoor  sound 
propagation,”  J.  Acoust.  Soc.  Am.  94, 2343-2352  (1993). 

Xiao  Di  and  Kenneth  E.  Gilbert,  "An  exact  Laplace  transform  formulation  for  a  point 
source  above  a  ground  surface,"  J.  Acoust  Soc.  Am..  93,  714-720  (1993). 


2.  Papers  submitted  to  refereed  journals 

K.  Attenborough,  S.  Taherzadeh,  H.E.  Bass,  R.  Raspet  G.R.  Becker,  A.  Gudesen.  A. 
Chrestman.  G.A.  Daigle.  A.L'sperance.  Y.  Gabillet  K.E.  Gilbert,  Y.L.  Li,  M.J. 
White,  P.  Naz.  J.M.  Noble.  H,  van  Hoof,  "Benchmark  cases  for  outdoor  sound 
propagadon  models,"  submitted  to  the  Journal  of  the  Acousdcal  Society  of  America 
(This  paper  involves  an  intemadonal  collaboradon  on  an  extensive  set  of  benchmark 
calculations.  The  names  are  in  alphabetical  order.) 


3.  Books  or  chapters,  authored  or  edited 

Kenneth  E.  Gilbert  and  Lucy  J.  Ameling,  "Surface  scattering:  the  role  of  turbulent  mixing 
of  microbubbles,"  to  be  publish^  in  a  refereed  hard-cover  book  of  the  Proceedings 
at  the  International  Conference  on  Theoretical  and  Computational  Acoustics,  Mystic 
Connecticut,  July  5-9,  1993. 

4.  Technical  reports  written 

Kenneth  E.  Gilbert  and  Ralph  R.  Goodman.  ".Acoustic  backscatter  from  the  near-surface 
oceanic  bubble  layer,"  contribution  to  Applied  Research  Laboratory  Review  (1993). 

Xiao  Di  and  Kenneth  E.  Gilbert  ".An  analytic  conformal  mapping  for  sound  propagation 
over  terrain,"  NCPA  Repon  No.  XD  1-93,  University  of  Mississippi  (Jan.  1993) 


5.  Invited  papers  or  talks 

Xiao  Di  and  Kenneth  E.  Gilbert.  ".Application  of  numerical  grid  generation  to  sound 
propagation  over  irregular  terrain."  invited  paper  present^  at  the  124th 
meeting  of  the  Acoustical  Society  of  America.  New  Orleans.  Louisiana.  3 1  October 
-  4  November,  1992. 

6.  Contributed  papers  or  talks 

Kenneth  E.  Gilbert  "Stochastic  scanering  model  predictions  for  reverberation  from  the 
near-surface  oceanic  bubble  layer."  presented  at  the  124th  meeting  of  the 
.Acoustical  Society  of  Amenca.  .New  Orleans.  Louisiana.  31  October  -  4  November, 
1992. 
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Xiao  Di  and  Kenneth  E.  Gilbert  "A  phase  screen  approach  to  sound  propagation  through 
small-scale  turbulence,  presented  at  the  124th  meeting  of  the  Acoustical  Society  of 
America,  New  Orleans,  Louisiana,  31  October  -  4  November,  1992. 

Kenneth  E.  Gilbert  and  Xiao  Di,  "Some  recent  developments  in  parabolic  equation  methods 
for  outdoor  sound  propagation,"  presented  at  the  26th  panel  meeting  of  the  NATO 
Research  Study  Group  11,  Quebec.  Canada.  October  26-29, 1992. 

Kenneth  E.  Gilbert  and  Lucy  J.  Ameling,  "Surface  scattering:  the  role  of  turbulent  mixing 
of  microbubbles,"  presented  at  the  Intemadonal  Conference  on  Theoretical  and 
Computational  Acoustics,  Mystic  Connecticut.  July  5-9, 1993. 

Xiao  Di,  Kenneth  E.  Gilbert,  and  Lucy  J.  Ameling,  "A  fast  phase-screen  method  for  sound 
propagation  through  a  turbulent  atmosphere,"  125th  meeting  of  the  Acoustical 
Society  of  America,  Ottawa.  Canada.  17-21  May  1993. 

7.  Conferences  attended 

124th  meeting  of  the  Acoustical  Society  of  America,  New  Orleans.  Louisiana,  31 
October  -  4  November,  1992. 

NATO  Research  Study  Group  1 1, 26th  panel  meeting,  Quebec,  Canada.  October 
26-29.  1992. 

Benthic  Boundary  Layer  Special  Research  Program  Workshop  I.  NRL-Siennis 
Space  Center,  MS.  October  29-30,  1992.  Sponsored  by  Office  of  Naval  Research 

Benthic  Boundary  Layer  Special  Research  Program  Workshop  n.  NRL-Stennis 
Space  Center,  MS,  January  7-8. 1993.  Sponsored  by  Office  of  Naval  Research 

Critical  Sea  Test  7  Data  Review  Workshop,  Applied  Physics  Laboratory.  Johns 
Hopkins  University,  February  3, 1993.  Sponsored  by  Office  of  Naval  Research 

125th  meeting  of  the  Acoustical  Societv  of  America.  Ottawa.  Canada.  17-21  May 
1993. 

Review  of  Army  Battlefield  Acoustics  Program.  May  24-24, 1993  Penn  State 
University.  Sponsored  by  the  Armv  Research  Laboimoiy  and  the  Army  Research 
Office. 

8.  Patent  applications 

None 

9.  Honors  (promotions,  awards,  etc.) 

Elected  Chairman  of  the  Technical  Committee  on  Underwater  Acoustics 

Elected  Vice  President  of  the  Pennsylvarua  Chapter  of  the  Acoustical  Society  of 
.America 

Appointed  to  membership  on  the  Technical  Committee  on  Acoustical 
Oceanography. 
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Organizer  and  session  chairman  tor  two  special  sessions  of  the  126th  meeting  of  the 
Acoustical  Society  of  America 

10.  Graduate  Students 

Name:  Lucy  J.  Kulbago 
Citizenship:  U.S. 

Tvpe  of  thesis:  theoretical 

Thesis  tide:  Acoustic  backscatter  hrom  the  near-surface  oceanic  bubble  layer 
Objective:  Incorporate  into  an  existing  plane-wave  surface  scattering  model  the 
improvements  necessary  to  model  backscatter  in  high  sea  states  and  shallow  water 
where  the  insonifying  field  departs  significantly  from  a  plane  wave.  Work  to  date 
has  dealt  with  high  sea  states  since  experimental  data  is  available. 

Name:  Timodiy  J.  Kulbago 
Cidzenship:  U.S. 

Type  of  thesis:  theoretical/numetical 

Thesis  title:  Wavelet  analysis  of  acoustic  backscatter  from  the  benthic  boundary 
layer. 

Objective:  Use  wavelet  analysis  of  acoustic  backscatter  from  shallow  water  bottoms 
to  extract  additional  information  relative  to  conventional  processing.  In  particular 
investigate  use  of  frequency  dispersion  characteristics  to  identify  s^ment  types 
and  to  identify  manmade  scatterers  in  the  bottom. 

Name:  Rodney  R.  Kone 
Citizenship:  U.S. 

Type  of  thesis:  theoretical/numetical 

Thesis  title;  A  fast  Green's  function  algorithm  for  sound  propagation  in  range- 
dependent  ocean  environments 

Objective:  Generalize  a  Green's  function  method  developed  originally  for 
atmospheric  acoustics  (propagation  in  a  half-space)  for  application  to  ocean 
acoustics  (waveguide  propag^on)  with  the  inclusion  of  the  effects  of  elasticity  and 
porosity.  Because  of  the  its  speed.the  method  should  be  extremely  useful  for  high- 
frequency,  broadband  calculations  required  in  shallow  water  acoustics  and  mine 
countermeasures. 

Name:  Scon  D.  Hansen 
Citizenship:  U.S. 

T>pe  of  thesis:  experimental 

Thesis  title:  An  atmospheric  turbulence  model  derived  from  inversion  of  sound 
pn^gation  measurements  in  an  acoustic  shadow  zone. 

Objective:  Make  accurate  measurements  of  the  sound  field  in  an  acoustic  shadow 
zone  where  the  sound  field  is  due  primarily  to  scattering  from  atmospheric 
turbulence.  Use  generalized  non-linear  inversion  techniques  to  extraa  an 
acoustically-bas^  empirical  model  for  near-ground  turbulence  in  the  atmosphere. 


11.  Postgraduate  Students 

None 

12.  Transitions 
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The  oceanic  bubble  layer  scattering  model  developed  by  Dr.  K.E.  Gilbert  on  6. 1 
binding  in  FY92  and  FY93  has  been  incorporated  by  Dr.  D.F.  McCammon  into  the 
operational  reverberation  model  used  by  the  Naval  Air  Weapons  Center  (  Hanh 
Environment  Low  Frequency  Reverberation  Model).  As  further  advances  are  made 
on  the  6. 1  research  model,  the  improven^nts  will  be  incorporated  into  the 
operational  model. 


BROADBAND  OCEAN  ACOUSTIC  TOMOGRAPHY 
(N00014-90-J-1365) 

Leon  H.  Sibul 

The  Pennsylvania  State  University 
Applied  Research  Laboratory 
P.O.  30,  State  College,  PA  16804 
(814)863-0312 
LHS2(8PSU.EDU 


OBJECTIVE 

To  investigate  the  use  of  the  broadband  ambiguity  function,  the  wavelet  transform  and 
tomx^raphic  concepts  for  characterization  of  time-varying  ocean  m*rdUrm  inhomogeneities. 

To  extend  classical  acoustic  tomographic  concepts,  that  essentially  assume  a  time  invariant  or 
a  slowly  varying  medium,  to  tracking  the  location  and  velocities  of  ocean  inhomogeneities. 
The  technique  uses  backscattering  of  broadband  active  acoustic  signals  from  moving 
inhomogeneities.  The  resulting  inverse  problem  is  solved  by  inverting  multidimensional 
wavelet  transforms. 

PROGRESS: 

We  have  shown  that  many  important  imaging  problems  including  radar,  sonar,  medical 
imaging  of  distributed  objects,  remote  srasing  of  the  earth,  oceans,  and  atmosphere, 
tooK^raphy,  and  time-varying  channel  identifictuion  fall  into  the  conceptual  firework  of 
the  theory  of  generalized  inverses  of  linear  operators.  The  general  problem  can  be  stated  in 
terms  of  the  following  operator  equation 

7k  =  r,  (1) 

where  r  is  the  received  or  sensed  signal  and  T  is  a  linear  operator  that  acts  on  x  (an  unknown 
image  or  medium  characterization).  T  is  described  by  a  known  measurement  system  or 
inobing  signal.  The  problem  is  to  invert  this  operator  equation  to  determine  the  unknown 
medium  characterizatitHi  x.  However,  the  inverse  of  r  may  only  exist  in  the  sense  of  a 
generalized  invme.  We  show  that  a  fundamental  theorem  for  generalized  inverses 
establishes  the  basic  processing  structure  for  remote  sensing  and  tomogrty^y  prd)lems. 

While  keqnng  the  generality  of  inverse  problems  in  mind  we  have  focus^  our  research  on 
estimation  of  broadband  and  narrowband  (NB)  spreading  functions  for  characterization  of 
randomly  time-varying  propagation  and  scattering  media.  The  key  step  in  processing  is 
calculation  of  the  adjoint  of  the  operator  T. 

The  adjoint  operator  associated  with  the  NB  delay-Doppler  spreading  function  estimation  is  a 
generalized  Gabor  trantform  of  the  received  signal  with  respect  to  the  transmitted  signal.  In 
the  widdtaml  case,  the  adjoint  operator  is  a  wavelet  transform  of  the  received  signal  with 
respe^  to  tlw  transmitted  signal  (mother  wavelet).  In  both  cases,  the  processor  output  is  a 
g^p  theoretic  convolution  of  the  spreading  function  with  the  appropriate  autoambiguity 


tWiMM 


Ic. 

Figure  1.  Imaging  and  Parameter  Estimation  of  a  Rough  Rotating  Sphere. 


function.  We  used  the  resolution  of  identity  (or  Shur-Frobenius-Godement  theorem)  to  show 
that  autoambiguity  functions  are  reproducing  kernels  for  some  subspaces  of  finite  energy 
sqneading  functions.  This  fact  is  used  to  design  probing  signals  that  can  give  good  estimates 
of  apprqniately  limited  spreading  functions  without  forming  the  inverse  operator.  Required 
signid  resolution  properties  are  determined  from  above  fundamental  theoretical  results. 

The  above  described  theoreticai  results  were  ai^iied  to  imaging  and  parameter  estimation  of 
rough  rotating  spheres  and  layers  of  randomly  moving  distributed  scatterers.  Rough  rotating 
^)heres  afproximaJbe  many  moving  inhomogendties  such  as  rotating  eddies,  clouds  of  bubbles 
and  other  objects  in  the  medium  that  are  observed  from  a  moving  platform.  Key  processing 
stqn  are  illustrated  on  Figure  1.  An  overhead  view  of  the  wideband  spreading  function  of  a 
rough  rotating  s{diere  is  shown  on  Figure  lb;  it  is  imaged  by  a  wideband  V-FM  signal  that 
has  the  autowavelet  transform  (autoambiguity  function)  shown  on  Figure  la.  The  observed 
image  at  die  0  dB  SNR  is  shown  on  Figure  Ic.  We  have  developed  a  generalized  Hough 
transform  technique  to  accurately  estimate  the  parameters  of  the  rotating  sphere  from  the 
noisy  and  blurred  image  (Figure  Ic). 
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PAPERS  PUBUSHED  IN  REFEREED  JOURNALS 

[1]  J.  A.  Tague  and  L.  H.  Sibul,  "Matrix  Representation  of  Propagation  and  Scattering 
Opmtions,"  Signal  Processing,  Vol.  29,  No.  3,  December  1992,  pp  299>307. 

[2]  R.  L.  Culver,  L.  H.  Sibul,  D.  W.  Ricker,  D.  F.  McCammon,  "High  Frequency 
Requirements  for  Shallow  Water  ASW,"  U.S.  Navy  Journal  of  UnderMOer  Acoustics, 
October  1992. 

[3]  M.  P.  Mahon,  L.  H.  Sibul,  and  H.  M.  Valenzuela,  "A  Sliding  Window  Update  for  the  Basis 
Matrix  of  die  QR  Decomposition,"  IEEE  Ttansaaions  on  Signal  Processing,  Vol.  41,  No.  S, 
pp  1951-1953,  May  1993. 

[4]  R.  K.  Young,  L.  G.  Weiss,  and  L.  H.  Sibul,  "Multisensor  Processing  for  Underwater 
Vehicles,"  /.  Underwater  Acoustics,  Vol.  43,  No.  1,  pp  167-178,  January  1993. 


accepted  for  PUBUCATION  in  REFEREED  JOURNALS 

{!]  A.  N.  Mirkin  and  L.  H.  Sibul,  "Application  of  die  Maximum  Likelihood  Criterion  to 

Multiple  Source  Localization  in  an  Inhomogeneous  Medium,"  acc^ted  for  publication  in 
J.  Acoust.  Soc,  Am.,  Fdiniary  1994. 

[2]  L.  H.  Sibul,  L.  G.  Weiss,  and  T.  L.  Dixon,  "Characterization  of  Stochastic  Propagation  and 
Scattering  via  Gabor  and  Wavelet  Transforms,"  accepted  for  publication  in  the  Journal  of 
Computational  Acoustics,  June  1994. 

[3]  L.  G.  Weiss,  R.  K.  Young,  and  L.  H.  Sibul,  "Wideband  Processing  of  Acoustic  Signals 
Using  Wavelet  Transforms,  Part  I;  Theory,"  to  appear  in  J.  Acoust.  Soc.  Amer. 


CHAPTERS  IN  BOOKS 

[1]  N.  K.  Bose  and  L.  H.  Sibul  "Multidimensional  Signal  Processing:  Sensor  Array  Processing," 
cli^iter  in  Electrical  Engineering  Handbook,  CRC  Press,  Inc.,  1993. 

[2]  L.  H.  Sibul,  R.  K.  Young  and  M.  L.  Fowler,  "Optimum  Signal  Processing  and  Medium 
Characterization  in  the  Time-Frequency  and  Wavdet  Transform  Domains,"  in  Acoustics 
Signal  Processing  for  Ocean  Exploration,  edited  by  Josd  M.  F.  Moura  and 

Isabd  M.  G.  Lourtie,  Kluwer  Academic  Publishers,  Dordrecht,  Netherlands  1993. 
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INVITED  TALKS 


[1]  L.  H.  Sibui,  ”Wavel«t  Transforms  and  Their  Applications,”  University  of  Tartu,  February  11, 
1993,  Taitu,  Estonia. 

[2]  L.  H.  Sibui,  "Wavelet  Transforms  and  Their  Applications,”  Technical  University  of  Tallinn, 
Fd)niary  17,  1993,  Tallinn,  Estonia. 

[3]  L.  H.  Sibui,  "Signal  Processing  -  Detection,  Tliresholding,  Signal  Ambiguity  Functions, 
Adaptive  Algoriduns,  and  Parameter  Estimation,”  presented  at  Advanced  Defense 
Development  Laboratory,  May  10-15,  1993,  Chinhea,  Korea. 


CONTRIBUTED  TALKS  (Also  published  in  Conference  Proceedings) 

[1]  T.  L.  Dixon*,  L.  H.  Sibui,  and  R.  K.  Young,  "Systems  Identification  via  Regularization  of 
the  Wideband  Inverse  Problem,”  Proceedings  of  Conference  on  Ii^rmation  Sciences  and 
Systems,  Johns  Hopkins  University,  Baltimore,  MD,  March  24,  1993. 

[2]  A.  P.  Chaiyasena*  and  L.  H.  Sibui,  "Group  Representation  Theory  and  Signal  Ambiguity 
Functions,"  Proceedings  of  Conference  on  Iiformation  Sciences  and  Systems,  Johns  Hopkins 
University,  Baltimore,  MD,  March  24-26,  1993. 

*  Presenter 


CONFERENCES  ATTENDED 

[1]  T.  L.  Dixon,  L.  H.  Sibui  and  A.  P.  Chaiyasena  atteaded  Conference  on  Information  Sciences 
and  Systems,  Johns  Hopkins  University,  Baltimore,  MD,  March  24-26,  1993. 
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GRADUATE  STUDENTS 


1.  A.  Peter  Chaiyasena . 

Ph.D  io  Matthonatics 

"Group  Theoretic  Foundations  of  Wavelet  Transforms  and  Signal  Ambiguity  Functions" 
OfeWfliYCi 

Devel<^  unified  group  theoretic  foundation  to  wavelet  and  other  group  dieoretic  transforms,  wide- 
and  narfDW-4)and  ambiguity  functions  and  Wigner  transforms.  An  important  conc^t  that  we  exploit 
in  this  study  is  diat  both  the  narrowband  and  widd)and  ambiguity  fon^ons  are  coefficients  of  the 
unitary  representttions  of  their  respective  groups.  Wide-band  ambiguity  functions  are  coefficients  of 
the  af^  group  and  narrowband  ambiguity  functions  are  coefficients  of  the  Heisenberg  group.  This 
foct  provides  insight  into  inqtortaitt  concepts  of  admissible  signals,  ambiguity  cons^ation,  and 
ambiguity  function  invariance  properties  foat  are  important  for  signal  design.  This  study  uses 
contractions  of  the  groups  to  study  the  relations  between  Heisenberg  and  affine  groups  from  which 
s^^ximations  of  wideband  ambiguity  functions  by  computationally  efficient  narrowband  ambiguity 
functions  can  be  established.  Theoretical. 

Student  Nationality:  American 
Graduated,  May  1993 


2.  Kenneth  L.  Hillsley 

Ph.D  in  Electrical  Engineering 

"Wavelet  Transform-based  Widd>and  Tomographic  Reconstruction  of  Stochastic  Media" 
Objective: 

The  objective  of  this  research  is  to  solve  the  inverse  problem  of  reconstructing  position  and  velocity 
ctmiptmeotB  of  turbulmit  inhomogeneities  by  developing  a  widd>and,  wavelet  transform-based, 
tomograidiic  reconstruction  theory.  By  considering  widd>and  processing  with  the  pownfol  technique 
of  tomogrqriiic  imaging,  it  will  be  possible  to  obtain  high-resolution  reconstruaions  of  multiscale 
time-v«ying  inhomi^eneities.  The  imerplay  between  wideband  processing  and  the  wavelet  transform 
is  an  emergiiqt  area  of  research.  The  unique  ability  of  wavelet  transforms  to  efficiently  and 
accurately  diaracterize  widdumd  and  nonstationary  signals  over  broad  parameter  ranges  make  wavelet 
transforms  a  naniral  analytical  tool  for  the  reconstruction  problem.  Theoretical  and  numerical. 

Student  Nationality:  Anwrican 
Expected  date  of  graduation:  May  199S 
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3. 


Tereu  L.  Dixon 
Ph.D  in  Electrical  Engineering 

"Wideband  Imaging  of  Distributed  Objects  Using  Wavelet  Transforms  and  Generalized 
Inverse  Theory" 

Obisstbo: 

The  objective  of  this  research  is  to  use  the  theory  of  generalized  inverses  and  wavelet  transform 
dieory  to  establish  optimum  processing  structure  for  active  acoustic  remote  imaging  of  time  varying 
distrUMited  objects,  and  randomly  time-varying  propagation  and  scattering  channels.  Results  from 
wideband  ambiguity  function  and  wavelet  transform  theory  are  used  to  establish  theoretical  limitations 
on  imaging  and  paramettf  estimation  of  distributed  objects.  Theoretical  results  are  verified  by 
computer  esperiments.  Theoretical  and  numerical. 

Student  Nttionality;  American 
Expected  date  of  graduation:  August  1994. 


TRANSITIONS 

Results  of  diis  research  are  continuously  transitioned  to  the  Shallow  Water  Torpedo  Detection  and 
Qassification  6.2  Program  and  the  Shallow  Water  Advanced  Technology  Demonstration  (6.3A) 
Program.  An  estimator/correlator  signal  processing  technique,  that  was  initiaiiy  researched  in  this 
program,  will  be  implementeu  and  tested  in  the  Shallow  Water  Detection  and  Classification  project. 
Reseirdi  results  from  this  project  are  used  for  die  Mine  Countermeasure  Program,  condition  based 
maintenance  and  ntm-destructive  testing.  Leon  Sibul  has  conducted  many  seminars  and  workshops, 
and  given  presentations  to  transition  research  results  for  the  6.2  community. 

Former  students  who  are  applying  resell  results  to  industrial  and  Navy  projects  are: 


Dr.  Mark  L.  Fowler 
Dr.  Adam  N.  Miridn 
Dr.  Dennis  W.  Richardson 
Dr.  }(dm  A.  Tague 
Dr.  Guy  R.  Sohie 


IBM  Federal  Systems 

NUWC,  Newport,  RI 

IBM  Federal  Systems 

University  of  Ohio,  Athens,  OH 

GE  Corporate  Research  Center,  Schenectady,  NY 
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A  Mid-Frequency  Sound  Speed  Measurement  Instrument 
N00014-9a-J-1365 

Ralph  R.  Goodman,  William  Thompson,  Jr. 

Applied  Research  Laboratory 
P.  O.  Box  30 
State  College,  PA  16804 
Telephone:  (814)  863>8140,  (814)  863-41S6 


An  iq^noach  to  the  problem  of  real-time  monitoring,  with  fine  spatial  resolution,  small  changes 
in  the  speed  of  sound  in  water,  such  as  might  be  associated  with  air-bubble  filled  water  near  the 
ocean’s  surface,  is  being  investigated.  The  thought  is  that  the  input  electrical  impedance  or 
admittance  of  an  electroacoustic  transducer,  which  is  a  function  of  the  acoustic  radiation  loading 
on  its  active  face,  can  be  made  an  even  more  soisitive  fimction  of  that  loading  by,  in  effect, 
amplifying  the  loading.  This  can  be  accomplished  by  locating  the  transducer  at  one  end  of  a 
planar  (one  dimensional)  wav^uide,  filled  with  the  water  to  be  evaluated,  and  terminated  at  the 
othn  end  by  some  impedance.  The  input  impedance  of  the  wav^uide  then  constitutes  the 
radiatum  load  on  the  transducer. 

Initially  it  was  thought  that  a  one-quarter  wavelength  long  waveguide  terminated  with  a  very 
large  impedance  would  provide  the  most  sensitive  situation  for  observing  small  changes  in  the 
sound  ^eed  of  the  ottrained  water.  However,  some  experiments  using  a  small,  but  low 
frequency,  beiuter-type  piezoceramic  transducer  positioned  at  one  end  of  a  thick-walled, 
cylindiick,  aluminum  tabe  which  was  terminated  at  the  end  opposite  the  transducer  by  a  thick 
aluminum  plate,  once  again  revealed  the  impracticality  of  trying  to  create  a  rigid  boundary  in 
a  fluid  me^um.  An  unusual  observation  during  these  experiments  was  the  tqq)earance  of  a 
second  resonance,  i.e.,  a  second  peak  in  the  input  conductance,  at  a  frequency  lower  than  the 
fundamental  mechanical  resonance  of  the  transducer  without  the  waveguide.  Since  this 
addititmal  resonance  is  wholly  an  artifact  of  the  presence  of  the  waveguide,  the  possibility  of 
expknting  it  as  a  barometer  of  changes  in  the  water’s  sound  speed  is  being  examined. 

A  madumnatical  model  of  the  wav^uide-loaded  transducer  has  been  developed  and  calculations 
of  the  input  conductance  of  the  transducer  do  reveal  two  resonances,  the  lower  frequency  one 
being  a  consequence  of  die  mass-like  reactance  of  the  waveguide  input  impedance  (assuming  the 
termination  impedance  of  the  waveguide  is  also  mass-like,  as  probably  was  the  case  for  the 
experintents  alluded  to,  and  assuming  the  loigth  of  the  waveguide  is  a  small  fraction  of  a 
wavelength,  as  was  the  case)  canceling  the  sdffrtess-like  impedance  of  the  transducer  at  a 
frequency  bdow  the  fiindamaital  mechanical  resonance.  More  importantly,  further  calculations 
show  that  the  frequoicy  of  this  auxiliary  resonance,  a  quantity  that  should  be  able  to  be 
measured  quite  rea^y,  varies  as  much  as  19%  for  a  20%  variation  in  the  ^)eed  of  sound  of  the 
entrained  water  and  in  a  very  linear  manner. 

The  cylindrical  wav^uide  used  in  the  experimoits  actually  had  a  number  of  axially  orirated 
^as  cut  dum^  the  wall  of  the  tube  to  permit  free  flooding  of  the  tube.  These  slots  seemed 
to  have  no  consequence  upon  the  existence  of  this  additional  resonance  but  would,  of  course, 
permit  the  wmer  to  change  as  an  apparatus  was  towed  through  the  ocean.  The  largest  of  the 
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small  bender  type  transducers  used  had  a  fundamental  resonance  at  about  2500  Hz  and  the 
auxiliary  resonance  occurred  as  low  as  1400  Hz.  The  entire  transducer  and  waveguide  assembly 
fit  in  a  cylindrical  ravelope  about  three  inches  in  diameter  by  no  more  than  five  inches  in 
length. 
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Problem 

The  propagation  of  acoustic  waves  through  turbulent  fluid  media  has  attracted  a  good 
deal  of  theoretical  attention.  This  is  so  because  of  the  number  of  areas  of  engineering 
interest  where  this  phenomena  is  present.  As  an  example,  in  undersea  applications,  the 
distorticHi  of  acoustic  wave  fronts  via  turbulence  has  a  clear  relevance  to  the  operation  of 
sonar  arrays.  The  processes  of  refraction,  advection,  and  scattering  in  a  moving,  turbulent 
flow  field  can  groszy  distort,  even  render  completely  incoherent,  an  incident  acoustic  wave. 
As  such,  an  understanding  of  the  mechanisms  of  this  distortion  in  the  context  of  turbulence  is 
required.  Early  cm,  the  interest  in  such  phenomena  was  primarily  concerned  with  the  effects 
of  laige  scale  turbulence  on  atmospheric  and  oceanic  acoustics.  In  terms  of  theoretical 
treatment,  most  of  the  early  work  involved  a  statistical  approach  to  the  problem  in  which  the 
effect  of  turbulence  was  accounted  for  by  a  modified  refractive  index.  More  recently,  a 
deterministic  approach  has  been  utilized  in  which  acoustic  waves  are  propagated  through 
instantaneous  realizations  of  turbulent  flows.  Statistics  are  then  determined  by  compiling 
results  from  a  lai:ge  number  of  such  realizations.  This  more  recent  work,  however,  has  been 
largely  limited  to  isotropic  turbulence.  Clearly,  for  practical  reasons,  it  would  be  useful  to 
extend  this  methodology  to  non-isotropic  turbulence.  That  is  the  ultimate  goal  of  this 
project. 

Approach 

In  this  woiic,  the  deterministic  approach  is  followed.  As  such,  the  problem  essentially 
becomes  a  two-part  problem:  first,  the  instantaneous  turbulent  flow  fields  are  generated,  and, 
sectmd,  the  acoustic  waves  are  computationally  propagated  through  a  number  of  these 
"hozen”  flow  fields.  Until  recently,  work  on  this  problem  has  concentrated  on  the  first  part. 
Specifically,  the  c^xdjility  has  been  developed  to  generate  two-  and  three-dimensional 
mlizations  of  homogeneous  isotropic  turbulence.  The  generation  process  essentially  begins 
with  die  specificatitm  of  a  turbulent  kinetic  energy  spectrum.  This  spectrum  is  then 
discretized  into  a  finite  number  of  wave  modes.  An  orientation  for  each  wave  vector  is  then 
randomly  assigned  according  to  probability  densities  that  will  yield  a  uniform  distribution  in 
wave  number  qiace.  Complex  amplitudes  for  each  of  the  wave  vectors  in  the  discrete 
plectrum  are  then  randomly  chosen  such  that  the  incompressibility  and  reality  con.>traints  are 
satisfied.  Longitudinal  and  lateral  correlation  functions  compiled  from  25000  such 
realizations  are  shown  in  figure  1. 

Currently,  work  cm  the  project  is  centered  on  the  second  aspect  of  the  problem; 
namdy,  the  acoustic  wave  propagation.  The  equation  governing  the  propagation  of  acoustic 
distorbaitees  in  a  randomly  inhomogeneous  medium  is  the  so-called  stochastic  Helmholtz 
e^iatkm.  Two  common  jqqnoaches  to  the  solution  of  this  equation  ate:  the  geometric 
aoou^cs  sqiproximation  (allowing  for  the  use  of  ray  tracing  techniques)  and  the  parabolic 
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^roximation.  Both  approaches  are  currently  being  implemented.  Presently  then,  the  main 
thrust  is  the  verification  of  these  basic  elements  against  previous  work  (analytical, 
computational  and  experimental). 

Future  Work  (to  be  completed  with  current  funding) 

As  stated  in  the  first  paragraph,  the  primary  objective  of  this  work  is  the  extension  of 
these  deterministic  techniques  to  non-isotropic  turbulent  flows.  This  will  necessarily  require 
a  complete  rethinking  of  the  first  part  of  this  work,  namely,  the  generation  of  the  individual 
realizations  of  the  turbulence.  The  intention  is  to  utilize  the  so-called  maximum  entropy 
family  of  probability  density  functions  to  construct  random  velocity  fields  satisfying  an 
arbitrary  number  of  moments  of  the  velocity.  With  this  approach,  for  example,  given  a 
Reynolds  stress  tensor  for  an  arbitrary  turbulent  flow,  it  would  be  possible  to  construct 
individual  realizations  of  the  turbulence.  Upion  constructing  a  large  number  of  such 
realizations,  and  ensemble  averaging  the  moments,  the  specified  Reynolds  stress  tensor  will 
be  recovered.  In  this  way  it  is  felt  that  information  from  many  standard  turbulence  models 
could  be  used  for  construction  of  the  flow  field  realizations.  The  number  of  moments  to  be 
specified  will  then  be  related  to  the  level  of  closure  of  the  particular  turbulence  model. 


Longitudinal  it  Loterol  Correlotion  Functions 
25000  reoiizotions 
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Figure  1.  Longitudinal  and  lateral  correlation  functions  for  25000  realizations  of 
homogeneous,  isotropic  turbulence  computed  using  random  Fourier  modes. 
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Prediction  of  the  strength  (or  failure)  of  materials  has  always  been  one  of  the  most 
fundamental  elements  of  engineering.  A  vast  body  of  knowledge  exists  about  the  strength  of 
materials  as  well  as  the  failure  mechanisms  under  static  loads.  However,  structures  under 
relatively  low  static  loads  (compared  to  their  material  yield  strengths)  are  known  to  develop 
cracks  and  fail  after  a  period  of  time  when  structural  vibration  is  present.  Some  common 
examples  of  these  structures  are  gears,  airplane  engine  pylons.  Naval  aircraft  tailhooks  and 
landing  gear,  and  even  machine  tools.  In  recent  years,  the  use  of  ultrasonic  non-destructive 
testing  has  allowed  die  detection  of  small  cracks  in  materials.  But,  by  the  time  a  crack  is 
large  enough  the  reflect  an  ultrasonic  wave,  much  of  the  structural  strength  may  already  be 
lost.  This  research  develops  a  technique  for  monitoring  the  growth  of  cracks  in  a  material 
using  a  simple  but  precise  measurement  of  fatigue-related  changes  in  the  structiual  damping. 

Any  growth  in  the  number  of  dislocations  should  also  increase  the  mechanical 
damping  due  to  wave  scattering  by  the  dislocations  as  well  as  increase  the  electrical 
resistance.  Note  diat  at  absolute  zero  temperature  most  metals  can  form  a  perfect  crystal  and 
become  superconducting.  The  increases  in  electrical  and  mechanical  damping  as  well  as 
increases  in  material  stiffitess  should  be  usable  by  a  sensor  system  designed  to  monitor 
fatigue.  According  to  C.R.  Barrett,  W.D.  Nix,  and  A.S.  Tetelman,  The  Principles  of 
Engineering  Materiids,  (Prentice-Hall,  1973),  p258.,  the  movement  and  presence  of 
dislocations  both  generate  and  interfere  with  vibration  waves,  but  the  effect  of  the  drag  forces 
is  small  compared  to  other  mechanisms.  Therefore,  as  a  structure  in  the  process  of  failure  is 
subjected  to  vibrations,  die  time  to  failure  should  decrease  wdiile  die  mechanical  damping  and 
sdffiie&t  should  increase. 

A  simple  but  crude  mq?eriment  ^paratus  is  shown  in  Figure  1  below.  A  3/8"  by 
0.030"  by  30"  long  band  of  high  carbon  steel  rated  at  1125  lbs.  breaking  strength  is  notched 
leaving  only  about  1  mm.  One  end  is  rigidly  attached  while  the  other  end  is  placed  in  tension 
through  a  stiff  spring  in  die  igiparatus.  The  purpose  of  die  spring  is  to  keep  die  band  in  near 
constant  tension  as  die  band  stretches  in  the  notched  area.  Widi  no  vibration  applied,  the 
band  gmerally  lasts  many  hours  before  breaking  -  if  it  breaks  at  all.  With  a  small  amount  of 
random  vibration  qiplied,  bands  widi  identical  notches  and  tension  consistendy  break  between 
30  minutes  and  slighdy  over  2  hours  time.  Tension  is  held  nearly  constant  by  initially 
"tuning"  die  band's  fundamental  frequency  to  die  same  value  for  each  test  and  by  carefully 
producing  identical  notch  widths,  ^so  a  spring  can  be  seen  at  the  right  end  in  Figine  1 
udiich  helps  keep  toision  constant  during  die  fatigue  process. 


Apparatus  for  Measuring  Damping 


-Top  view  of  X)3(r  Band 


Shaker/Sllnoer  Accelerometer 

Figure  1  Damping  measurements  are  achieved  by 
observing  the  Q  of  structural  resonances 


A  broadband  random  noise  signal  is  applied  to  the  shaker  which  excites  the  band  via 
die  stinger  as  seen  in  Figure  1.  An  accelerometer  senses  the  response  of  the  band  during  die 
fatigue  test.  The  positions  of  the  shaker  and  accelerometer  have  a  profound  affect  on  the 
measured  transfer  function,  but  not  on  the  measured  changes  in  damping.  Figure  2  shows  two 
typical  transfer  function  frequency  responses  taken  just  after  initial  setup  and  after 
approximately  90  minutes  of  vibration  excitation. 


Figure  2  Very  small  increases  in  material  stiffness 
and  larger  increases  damping  are 
observable  over  the  90  minute  test. 
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eradicate  Students: 

1)  Mr.  Douglas  Koehn,  USA,  MS  graduation  expected  Fall  1994,  Experimental, 
"Vibration-Induced  Material  Fatigue  Measurement  Using  Modal  Damping" 

Objective:  The  thesis  defines  the  known  fundamental  physical  mechanisms  for  sonic 
fatigue  of  materials  and  their  relationship  to  damping.  A  simple  apparatus  is 
constructed  for  transfer  function  measurements  of  a  notched  band  under  static  and 
dynamic  loads.  Photo-microgri^ihs  compare  physical  observations  of  crack  formation 
and  growth  to  corresponding  structural  damping  measurements.  The  trend  in  damping 
changes  during  failure  will  be  evaluated  as  a  precursor  to  structural  failure. 


2)  Mr.  Sevag  Arzoumanian,  USA,  MS  graduation  expected  Fall  199S,  Experimental, 
"Active  Vibration-Induced  Fatigue  Mitigation" 

Objective:  This  thesis  iqiplies  active  vibration  control  using  adaptive  self-optimizing 
algoridims  to  stop  crack  growdi  in  a  seeded  fault  area.  Since  the  adaptive  control 
"error"  sensor  cannot  be  placed  directly  in  die  fault  area,  the  complete  dynamics  of  the 
notched  band  will  be  examined  to  develop  a  "remote  error"  sensor  system  consisting 
of  an  array  of  accelerometers.  The  response  of  the  adtqptive  controller  will  be  used  to 
simultaneously  monitor  structural  integrity  and  to  develop  a  material  condition 
prognosis  system  linked  widi  an  active  fatigue  control  mediodology. 
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